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ABSTRACT 

From new observations and literature data we investigate the presence of H I, dust, and optical cores in the central kiloparsec of 
low-power radio galaxies. The goal of this pilot study is to identify physical relations between these components, which can help 
us to study kinematics and feeding mechanisms in future samples of active galaxies. Our results are consistent with neutral gas 
being associated with dust on sub-kiloparsec scales. Objects that have HI absorption always have significant amounts of dust in their 
host galaxy. If there is no visible dust in the host galaxy, there is also no H I absorption. The presence of an unresolved optical core 
correlates with the H I column density, with the core being absent in high column density sources. This work opens a path for studying 
the kinematics of cold material in the central regions of active galaxies by combining information of H I absorption and molecular 
lines. Consistent with previous work, we find no evidence for a compact, parsec-scale obscuring torus in low-power radio galaxies. 
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1. Introduction 

The central regions of active galaxies (AGN) harbour essential 
clues to how these amazing objects work. Differences on small 
scales lead to enormous variety of observed characteristics on 
large scales. For radio-loud AGN, a distinction is made between 
high- and low -power radio galaxies, referred to as FR II and FR I, 
respectively (Fanaroff & Riley 1974) . Differences between these 
sources occur in radio morphology, optical, and X-ray proper- 
ties. They are attributed to different accretion rates or accretion 



modes (see e.g. Buttiglione et al. 2010 1. This implies differences 



in kinematics and morphology of gas and dust in the central re- 
gions. 

One difference between high- and low-power radio galax- 
ies is well established. Unified schemes postulate a parsec-scale, 
thick dust torus to explain the bimodality in optical spectra 
of high-power FRII radio galaxies (Barthel |1989 Antonucci 



1993| e.g.). In contrast, the optical nuclei of FRI radio galax- 



ies appear unobscured (e.g. de Koff et al. 1996 Chiaberge 



let al. 



have 



1999; |Capetti et al.|2002| i, though th e host galaxies o ften 



dloparsec-scale dust r i ngs and lanes (|de Koff et al.|2 000; 



Verdoes-Kleijn et al. 



2002| |de Ruiter et al.||2002| l. |ChTaberge 



et al. ( 1999 ) find that the brightness of the optical core correlates 



with the brightness of the radio core. Therefore, they conclude 
that the optical core emission is non-thermal, coming from the 
active nucleus. This is only possible if a parsec-scale torus is ab- 
sent. X-ray observations confirm this absence, and find evidence 
for a different accretion mode in FR I radio galaxies ( Balm averde] 
|et al.|2006||Hardcastle et al.T20Q9] l. 

A major question in AGN research is how the material is 
transported from kilo-parsec scales to the sub-parsec scale. This 
requires kinematical studies at high spatial and velocity resolu- 
tion. Due to their proximity, FR I radio galaxies are best suited 



for this purpose. HI absorption observations are a powerful tool 
for studying kinematics on small scales (e.g. van Gorkom et al. 
IT989}|Morganti et al.|200"T][20rj3||Vermeulen et al.|2003||Gupta 
|et al.|2006| l. Molecular lines can provide essential additional in- 
formation, as they also trace the cool interstellar medium. In this 
paper we present a pilot study to explore the relation between 
dust (associated with the molecular lines) and neutral gas (asso- 
ciated with the H I) on kilo-parsec scales or smaller. If we can 
relate the gas and dust on sub-kiloparsec scales, future studies 
can combine information from these sources to generate a com- 
plete kinematical picture. 

We expand upon earlier studies by combining H I informa- 
tion with HST imaging. Our primary goal is to determine if the 
presence of an optical core is associated with the absence of H I 
absorption. Next, we search for a relation between the dust and 
HI, where we expect to see absorption if dust lies in the line- 
of-sight to the radio source. If H I is indeed associated with the 
dust, deeper HI observations can be combined with molecular 
line data to study the kinematics of the central regions in a com- 
plete sample. This work is essential to provide new insight into 
AGN fueling and feedback. 



2. The sample 

Our sample consists of low-power radio galaxie^] logP\ 4ghz < 
25.5 W Hz . They all have typical FRI (core-dominated) radio 
morphologies. The host galaxies are early-type galaxies, with a 
range of ellipticities and environments. For our new VLA ob- 
servations we selected a representative sample of eight FR I ra- 
dio galaxies from the HST snapshot survey of the B2 catalogue 
( |Capetti et al.|2 002 ). The HST snapshot images trace dust down 



Throughout this paper we use Hq= 70 km s 1 Mpc 
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Table 1. VLA sample description (columns 1 through 6) and observational results (columns 7 through 10). 



Name 


ID 


z 


So,! A 


P (,4.9 


^■obs 


S abs 


rms 


Tim 


N m 


B1950 






(Jy) 


(WHz"l) 


(min) 


(mJy birr 1 ) 


(mJy birr 1 ) 


(kms-l) 


(xl0 20 cirT 2 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


0055+30 


NGC315 


0.0165 


0.464 


23.84 


27 


27.0 


1.0 


4.5±0.7 


8.3±0.3 


0104+32 


3C31 


0.0169 


0.067 


24.11 


66+47 


<1.8 


0.6 


<2.7 


<5.19 


0755+37 


3C189 


0.0428 


0.137 


24.62 


22 


<3.0 


1.0 


<2.2 


<4.23 


1322+36 


NGC5141 


0.0176 


0.054 


23.45 


57 


3.3 


0.6 


5.5+1.0 


10.0+1.8 


1346+26 


4C26.42 


0.0632 


0.135 


24.35 


139 


3.0 


0.5 


2.2+0.4 


4.1±0.7 


1350+31 


3C293 


0.045 


1.652 


24.93 


20 


148.0 


1.0 


12.8+0.1 


23.5±0.2 


1626+39 


3C338 


0.0303 


0.138 


23.98 


17 


<1.8 


0.6 


<1.3 


<2.5 


2116+26 


NGC7052 


0.0156 


0.040 


22.80 


88 


<1.2 


0.4 


<2.6 


<5.1* 



(2010) set N m < 0.34x 10 2 ' 



Notes. ( *' Upper limit from our VLA observations. Emonts et al. 
Clarification of the columns: (1) B2 rado name, (2) ID used in the paper, (3) redshift, (4) core flux at 1.4 GHz, (5) total power at 4.9 GHz, (6) total 
on-source integration time, (7) peak of the absorbed flux (8) rms noise, (9) integrated optical depth, and (10) HI column density. 
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Fig. 1. Plot of the continuum subtracted spectra for two of the 
detected sources from the VLA sample. The dashed vertical line 
indicates the systemic velocity. The x-axis is in optical barycen- 
tric velocity. The spectra are obtained against the brightest part 
of the radio source, which in both of these cases corresponds to 
the core. 



to scales of ~100pc and are used to identify the presence of an 
optical unresolved core. 

To increase the sample size, we have added 14 FRI radio 
galaxies from literature. They were selected to match our cri- 
terion for radio power, and to have HI data and optical HST 
images available. This literature sample consists of results ob- 
tained by |van Gorkom et al7] (| 1989| hereafte r JvG89), and four 
additional s ources from |van Gorkom et al. (19861^06 Young 
|etaLl(f T973] l, |Jaffe & McNamara| ( 11994] irand |van der Hulst et al. 



(1983 
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Fig. 2. Plot of the HI absorption detected in IC 1459 with the 
ATCA telescope. For comparison with the VLA spectra, the plot 
has been smoothed to a velocity resolution of 40 km s _1 . 



The total sample has 22 objects and is not meant to be com- 
plete, but representative of this type of radio sources. The full 
sample is listed in Table|2] 



3. Observations and H I results 

subsection VLA observations We observed eight sources with 
the VLA in A-configuration during December 2000 and January 
2001 . Observations were made with 6.25 MHz bandwidth and 64 
channels. The data were reduced with CASA^j We performed 
the default steps for HI absorption line data reduction, which 
are described in the online manual for H I absorption reduction 
in CASA0 Hanning smoothing was applied after the calibra- 
tion steps. Observational parameters for each source are listed 
in Table [1] 

H I absorption is detected in four out of eight sources (see 
Table [TJ. The peak absorbed flux exceeds 5cr, with cr being 
the rms noise per channel. The source 4C 26.42 is a new de- 
tection. The 3<x absorption peak in NGC 7052 is considered a 
non-detection, though it shows a hint of a resolved line. Since 
our observations, four sources have been studied with other tele- 
scopes. NGC 315 is presented in |Morganti"eTaL](p009); 3C 293 
was observed at high resolution by |Beswick et al.| "p002); and 
NGC 5 141 and NGC 7052 were observed at lower resolution in 



" http://casa.nrao.edu 

3 http://casaguides.nrao.edu 
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Table 2. H I and optical properties of the sample. First we list the sources with H I absorption, then we list the non-detections. The 
two rows labelled as 'Total' give the total number of sources that are listed above, and the number of objects that have either dust 
structures or an optical core. 



rS lyjv 


ID 




Ker 


Dust 


Diameter 


i 


r\ 
U 


optical 






1 f\20 „- m —2\ 

(X1U cm ) 




morph. 


(pc) 




( ) 


core 


(1) 


f~) \ 
W 




(4) 


(5) 




\i) 


w 


few 

(y) 


UUjjtjU 


IN V J v .1 i J 


a. j 


i 
i 


disk 




/ j 


n 1 s 


y 


0116+31 


4C31.04 


10.8 


2 


extended 


2586 


83 





n 


0238-08 


NGC1052 


3.0 


3 


disk 


416 


72 


0.2 


v 


0316+41 


3C84 


2.0 


4 


extended 


17625 


37 


1.0 


v 


1146+59 


NGC3894 


18.2 


2 


disk 


911 


74 


0.07 


n 


1216+06 


3C270 


7.2 


5 


disk 


263 


65 


0.07 


v 


1322-42 


Cen A 


50.0 


6 


extended 


3365 


90 







1322+36 


NGC5141 


10.0 


1 


disk 


638 


78 


o 


n 


1346+26 


4C26.42 


4.1 


1 


extended 


9063 


57 


0.14 


v 


1350+31 


3C293 


23.5 


1 


extended 


4130 


67 





n 


2254-367 


IC1459 


>1.2 


7 


disk 


310 


54 


o 


v 
j 


Total 




11 




11 








5 


UlUtTJZ 


i 




i 

i 


disk 


9 1 7ZL 


J 1 


W.o 


y 


0305+03 


3C78 


<1.47 


2 


none 








y 


0430+05 


3C120 


<1.21 


2 


none 








y 


0755+37 


3C189 


<4.2 


1 


none 








y 


1142+19 


3C264 


<1.65 


2 


disk 


675 


20 


0.2 


y 


1514+07 


3C317 


<1.70 


2 


extended 


2595 


33 


0.25 


y 


1626+39 


3C338 


<2.5 


1 


extended 


3463 




0.6 


y 


1637+82 


NGC6251 


<0.3 


1 


disk 


678 


65 


0.11 


y 


1652+39 


4C39.49 


<4.52 


2 


none 








y 


1807+69 


3C371 


<2.96 


2 


none 








y 


2116+26 


NGC7052 


<0.34 


8 


disk 


1227 


75 





y 


Total 




11 




6 








n 



Notes. Clarification of the columns: (1) B2 radio name, (2) ID used in this paper, (3) H I column density (4) reference for column 3, (5) morphology 
of the dust, (6) size of the dust structure, (7) inclination of the disk (where applicable), (8) distance between optical core and the dust detected in 
the optical images, and (9) presence of an optical core: y=yes, n=no. 



Re ferences. (1) This paper; (2 ) |van Gorkom et al.|jT989) ; (3) |van Gorkom et al.|(l985) ; (4) |De Young et al.| ( fT973> ; (5) Jaffe & M cNamara| ( [T994"l >; 
(6) |van der Hulst et al.|fl98ll >; (7) |Oosterloo e t al.| ( |1998| >; (8) |Emonts et al.|P010| >. 



Emonts et aT] ( |2010| >. Only NGC 7052 was not detected, and a 
stringent upper-limit was set. Our column densities and H I pro- 
files are consistent with literature results. With our limited band- 
width and sensitivity we cannot recover the very broad compo- 
nent in 3C293 (Morgant i et al.|2003| ). Two representative spec- 
tra, including 4C 26.42, are shown in Figure [T] 



3.1. Additional data 

A detection of HI absorption in IC 1459 has been obtained 
with the Australia Telescope Compact Array. IC 1459 was ob- 
served with 16.0 MHz bandwidth, at a velocity resolution of 
13 km s -1 (before Hanning smoothing), and with the 375 m con- 
figuration (Oosterloo et al. 1998 see Fig. [2]). The rms noise is 
1.4mJybeam _1 after Hanning smoothing. Although the spatial 
resolution of these observations is very poor, the background 
source is unresolved on arcsecond scales (Sad ler et al.| T989). 
Therefore, we can assume that the absorption is coming from 
these small scales. 



The giant radio galaxy NGC 625 1 is listed as a non-detection 
by JvG89. However, Evans et al. (2005 ) reported an unpublished 



column density of 4.5 X 10 /u cm _/ . To confirm this, we obtained 
a 12 hour WSRT service observation, with 1024 channels and a 
bandwidth of 20 MHz. The observation was processed similar to 
the VLA observations. We found no detection and set an upper 
limit of 7Y H i < 0.3 x 10 20 cirT 2 for this source. 



3.2. HI column density 

To derive the HI column density in units of cirT 2 , we used the 
standard equation 



JVhi = 1.835 x 10 1 



J r(v)dv • 



1.835 x 10 20 -t.s .AVoh,(l) 



in which T$ is the spin temperature and f c the covering fraction, 
t s is the optical depth summed over the channels where absorp- 
tion is present, and Av c h is the velocity width of a single channel. 
The product • Au c h is the integrated optical depth, r lnf . We as- 
sumed a spin temperature of 100K and a covering fraction of 
1. 

For non-detections we replaced the integral with the outcome 
for a Gaussian profile 



T(v)dv - 1 .06 ■ T p ■ AfpwHM 



(2) 



in which r p is the peak optical depth and Ai/fwhm the velocity 
width of the Gaussian at half maximum. In a non-detection the 
peak absorption depth is assumed to be less than 3<r, with cr 
equal to the rms noise per channel. Since r p * Sob s jS c , we can 
replace t p by 3<r/S c , with S c being the core flux at 1.4 GHz. The 
resulting upper limit of column density is then 



Nm < 1.93 x 10 



; 2() 



3cr 
5V 



(3) 
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no dust 
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HI detections HI non-detections 

(a) Histogram of the distribution of dust. 





HI detections HI non-detections 

(b) Histogram of the distribution of optical cores. 

Fig. 3. Histograms of the distributions of dust (a) and optical 
cores (b) in the H I detections and non-detections. 

We assumed the average velocity width of the detected 
sources for the line width Ai/fwhm = lOOkms -1 . For a proper 
comparison, we recalculated the literature upper limits, which 
use different values. 

The HI absorption technique is limited to objects with a 
bright core. Hence, in sources with a faint core, we may be 
missing HI if it has a low column density. In our small VLA 
sample this bias is not obvious. In the source with the weakest 
core, NGC 5141, H I is detected with an integrated optical depth 
of 5.5 + 1.0. The results for all VLA observations are listed in 
Table [1] 

Combining all the results, we have 1 1 detections and 1 1 non- 
detections. Owing to our selection method, these numbers can- 
not be used to derive detection rates for the population of low- 
power radio galaxies as a whole. 

4. Results 

The goal of this study is to compare the detection of H I in ab- 
sorption with the presence of kiloparsec-scale dust and an optical 
core. 

4.1. Relation between neutral gas, dust, and optical core 

The trends we find are illustrated in Figs. [3] and [4] and are listed 
in Table|2] They can be summarized as follows: 

1. All HI absorption detections have significant amounts of 
dust in their host galaxies, and sources without visible dust 
in the host galaxies have no H I absorption (Fig. 3(a) i. 

2. All sources without HI detection have an optical core 

(Fig.ro. 




20.0 20.5 21.0 

log N m (cm 2 ) 



22.0 



Fig. 4. Distribution of HI column density for the sample. The 
detections are grey, the detections with an optical core are 
hatched. Within the detections there is a clear trend for the 
sources without an optical core to have a higher H I column den- 
sity. This is reflected in the entire sample, taking into account 
that the non-detections in the white area only have upper limits 
of H I column density and all have an optical core. 

3. Sources with Nm > 10 21 cirr 2 have no optical core (Fig.jijl. 

We found sources with HI absorption that have an optical 
core, see Fig. [4] These sources have either extended dust far 
from the nucleus (see section 4.2), or disks with low inclination, 
consistent with a line of sight through a less dense medium. If 
Nh > 10 21 cirr 2 , the unresolved optical core is not detected. 
This column density corresponds to a visible extinction of ~ 1 
magnitude for Galactic extinction (Predehl & Sch mitt|[l995) . 
Because we have no handle on the intrinsic brightness of the 
optical core, future work on a larger, complete sample will have 
to clarify the scientific merit of our third result. 



4.2. Dust morphology 

In the optical images we find three distinct morphologies for the 
dust: 

- Extended: dust extending up to tens of kpc. In many sources 
this dust does not show any symmetry with respect to the 
core of the galaxy, often it is diffuse and/or filamentary. 
Some extended dust structures show evidence for smaller 
scale arms, or have a feather-like appearance. In a few 
sources it appears as a large dust lane (e.g. Centaurus A). 

- Disks: dust in a well-defined and symmetric distribution. The 
dust disks are rather compact, having radii of several hun- 
dreds of parsec. We find no disks with radii larger than 1 kpc. 

- None: in the HST images there is no visible dust. 

Our classification of dust structures is consistent with results 
from IdeRuiteretaL] pOOl) . 

The 1 1 sources with H I absorption all have dust, and are split 
between extended dust (5) and disks (6). The disks are at high 
inclinations, with a median value of 73°. Extended dust either 
covers the optical core, or is located at less than one HST reso- 
lution element from the core. The exception is 3C 84, where the 
dust is more than 1" from the optical core. This source is known 
to have a galaxy in the line of sight. This was first discussed by 
Min kowski! ( |1957[ ) and later confirmed by |Lynds| ( |1970| l. Two 
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separate H I absorption components are seen. The velocity dif- 
ferences are consistent with a component associated with the 
foreground galaxy and another with the radio galaxy (De Young 
|et al.||1973) |Rubin et aL|l977| ). It is not clear where the dust is 
located. VLBI observations are ongoing to spatially resolve the 
H I absorption components in this complex source. 

Among the 11 HI non-detections, 6 sources have dust. In 
four cases the dust is either in a face-on disk (e.g. 3C31) or 
at large projected distance from the core. The remaining two, 
NGC 6251 and NGC 7052, have an edge-on dust disk. The non- 
detection of NGC 625 1 can be explained by the existence of a 
strong warp in the inner regions. While the outer parts of the disk 
are nearly edge-on at 76° inclination, the inner regions are at 35° 
inclination ( Ferrarese & Ford|1999 1. If neutral gas is present in 
the dust disk, this would significantly suppress the column den- 
sity of H I in the center. NGC 7052 has the weakest radio core in 
our sample, which makes it the most difficult source to detect. 



Nevertheless, Emonts et al. ( 201Q) > set a stringent upper limit. 
Higher resolution observations of this source may provide in- 
sight into the structure of the disk. 

In sources without dust there is no detection of H I, and they 
have an unresolved optical core. Though this is consistent with 
the absence of obscuration on all scales, for these five sources 
alone we cannot fully exclude orientation effects. Considering 
the sample as a whole, we would expect a tighter relation be- 
tween the absence of an optical core and the presence of H I ab- 
sorption if a small-scale torus were present. This is not reflected 
in our results. 

5. Conclusions 

We have studied the presence of HI, dust, and optical cores 
in the central kiloparsec of low-power radio galaxies. We find 
that the presence of H I absorption is linked to the presence of 
dust. Furthermore, the absence of H I absorption is linked to the 
presence of an optical core. The optical core is not detected if 
A/hi > 10 21 crrT 2 . Our results are consistent with neutral gas 
being associated with extended dust and dust in the kiloparsec- 
scale disks. We require no parsec-scale dust torus. 

The distribution of H I may be more complicated than a reg- 
ular rotating central disk. This is evident from structure in the 
absorption profile, and is demonstrated by several higher resolu- 
tion studies. HI distributed in clouds with a filling factor much 
lower than unity, and H I in infalling and outflowing structures 
are seen (see the case of NGC315, Morgan ti et al.|20 09). Only 
VLBI observations will allow one to distinguish these effects. 
For sources with HI absorption, the line kinematics is an essen- 
tial component for understanding the feeding mechanism of ac- 
tive galaxies. Further insight into the physics of kiloparsec-scale 
disks will require deep observations of a complete sample. These 
are currently ongoing. More work will be done by the large sur- 
veys planned with the SKA pathfinders. 

Finally, molecular line transitions can provide crucial com- 
plementary information on the kinematics of the disks. This is 
especially important for sources that do not have a bright, unre- 
solved radio core. This work opens a potential path for combin- 
ing molecular lines with H I absorption to study kinematics on 
sub-kiloparsec scales. Systematic studies of the cold, molecular 
gas component are now starting. However, Flaquer et al. (2010) 
found that radio galaxies tend to have relatively low molecular 
gas masses. Given the limited sensitivity of current instruments, 
ALMA will make a major leap in this field. 
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